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Abstract 

The main aim of this paper is to study the scattering amplitudes in gauge field 
theories with maximal supersymmetry in dimensions D = 6,8 and 10. We perform a 
systematic study of the leading ultraviolet divergences using the spinor helicity and 
on-shell momentum superspace framework. In = 6 the first divergences start at 
3 loops and we calculate them up to 5 loops, in H = 8,10 the first divergences start 
at 1 loop and we calculate them up to 4 loops. The leading divergences in a given 
order are the polynomials of Mandelstam variables. To be on the safe side, we check 
our analytical calculations by numerical ones applying the alpha-representation and 
the dedicated routines. Then we derive an analog of the RG equations for the 
leading pole that allows us to get the recursive relations and construct the generating 
procedure to obtain the polynomials at any order of perturbation theory (PT). At 
last, we make an attempt to sum the PT series and derive the differential equation 
for the infinite sum. This equation possesses a fixed point which might be stable or 
unstable depending on the kinematics. Some consequences of these fixed points are 
discussed. 
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1 Introduction 


In the recent decade there was considerable progress in understanding the structure of 
the amplitudes (the S-matrix) in gauge theories in various dimensions (for review see, for 
example, m)- The gauge and gravity theories with maximal supersymmetries in D = 4, 
A/" = 4 SYM and A/" = 8 SUGRA are the most important examples. This progress became 
possible due to the development of the new techniques: the spinor helicity and momentum 
twistor formalisms, different sets of recurrence relations for the tree level amplitudes, the 
unitarity based methods for the loop amplitudes and various realizations of the on-shell 
superspace formalism for theories with supersymmetry [1]. 

The subject of investigation was mainly related to the so-called maximally supersym¬ 
metric theories, which are believed to possess special properties due the highest sym¬ 
metries. One of the insights was the discovery of the dual conformal symmetry for the 
AT = 4 SYM. Taking together the algebras of ordinary (super)conformal symmetry and 
the dual (super)conformal symmetry can be fused into an inhnite dimensional Yangian 
algebra [2] which in principle should completely dehne the S-matrix of the A/" = 4 SYM 
theory [3llll|5]. 

While the A/” = 4 SYM theory is completely on shell UV hnite and possesses only 
the IR divergences, in higher dimensions the situation is the opposite: there are no IR 
divergences even on shell but all theories are UV nonrenormalizable by power counting. 

It should be noted that the spinor helicity formalism and the unitarity based methods 
can be generalised to space-time dimension greater than D = The new 

computational methods gave new birth to the investigation of the UV properties of the 
S-matrices of formally nonrenormalizable gravity theories with extended supersymmetry 
(D = 4 A/" = 8 SUGRA is a particular example). The results obtained so far are in some 
sense controversial puniiii]. 

Among the gauge theories in higher dimensions with maximal supersymmetry there 
are the following four cases: 

D = 4Af = 4, D = 6Af = 2, D = 8 Af = 1, D = W Af = 1. 

No wonder if all these theories obey some exceptional properties. In this context, it is 
interesting to note that the integrands of the four-point amplitudes in any SYM theory 
have almost identical form (only the tree level amplitudes which are the common factors 
are different) and are heavily constrained by the dual conformal covariance in dimensions 
D < 6 [T^ (and likely in all dimensions T* < 10 [H]). 

The aim of this paper, which is a continuation of our previous papers HSHH] , is to 
investigate the amplitudes and their UV properties in maximally supersymmetric gauge 
theories in various dimensions. Namely, we evaluate the leading UV divergences in the 
four-point amplitude on shell in a number of loops and investigate their properties in all 
loops. The paper is organized as follows: in section 2, we briefly describe the spinor helicity 
formalism in H > 6, and in section 3, we consider the on shell superspace formalism in 
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D = 6,8. In section 4, we consider the structure of the colour ordered partial amplitudes. 
In section 5, we present the evaluation of the leading UV divergences both analytically 
and numerically. Section 6 summarizes the results of perturbative computation for various 
dimensions. Finally in section 7 we derive the all loop recursive relations for the leading 
divergences and make the attempts to summarize the whole PT series. The conclusion 
contains some speculations regarding the observed pattern and its implications on possible 
scenarios of UV hniteness of gauge and gravity theories. 


2 The spinor helicity formalism in various dimen¬ 
sions 

As was stated in the introduction, the spinor helicity and the on shell momentum su¬ 
perspace formalisms play a crucial role in the recent achievements in understanding the 
structure of the S-matrix of four dimensional supersymmetric gauge held theories. Here 
we discuss the generalization of these formalism to the case of even dimensions D = 6,8 
and 10. 

In our discussion we manly follow na. In even dimensions one can always choose the 
chiral representation of the gamma matrices as P^ (as usual {P^, P*^} = 27]^'^)-. 


PM _ 


0 




0 


( 1 ) 


where jj, is the SO{D — 1,1) vector representation index, A and B' = l,...,2^/^“^ are 
the Spin{SO{D — 1,1)) indices. We are interested in the cases D = 4,6,8,10. In this 
notation one can decompose the Dirac spinor -0 as a pair of Weyl chiral and anti-chiral 
spinors and Xa'- The Lorentz rotations of and A^' look like 


5Xa' = 


( 2 ) 


where 


cr 


liv\A _ 

)B — 




liu\B' _ 

)A' = 
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i 
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(a^)Am(cr") 


AB' 


CT 


(aY^ia^AB' 


(o'^)a'B 


(3) 


One can combine two Dirac spinors V ’2 and V’l into the Lorentz invariant combination 
using the charge conjugation matrix C dehned so that 

C'P'^C'-i = -(P^)^. 


(4) 

For the Weil spinors there are two possible decompositions of C depending on dimension: 

C = ( 1 for D = 4,8, (5) 
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and 

^=(f = (8) 

The r2 matrices obey the following relations 

VLba^^^ = 5%, VLb'A'^^'^' = SS' for D = 4,8, (7) 

and 

= 5%, = 6^1 for D = 6, 10. (8) 

The matrices hi can be used to raise and lower the indices of the spinors 

Aa = X^VtBA A^' = for Zl = 4, 8, (9) 

and to relate the chiral and antichiral spinors 

Aa = Aa', A^' = for D = 6, 10. (10) 

One can also construct the Lorentz invariants for the pair of spinors which are labeled by 
i and j: 

XfVtBA^f = (0), As/^iO^'^'AA'j = [ij] for D = 4,8, (11) 

and 

AA',iOA= W)^ XfVLfXA'j = {i\i] for D = Q, 10. (12) 

The matrices C can be always chosen in such a way that 

= -O for D = 4,10, 

= C'forT) = 6,8. (13) 

In some dimensions one can also construct additional Lorentz invariants. For example, 
in Z) = 6 one has Spin{SO{b,l)) = SU{4)*, so one can construct the invariants as 
contractions of spinorial indices A with the absolutely antisymmetric tensor Sabce asso- 
siated with SU{4)*-. eABcnXfX^X'^Xf = (1234), e^^^^XAA^B, 2 X 0 , 3 X 0,4 = [1234]. 

To relate the light like (massless) momentum with the pair of Weyl spinors we 
consider the Dirac equations for the spinors A"^ and Aa': 

(p^a^)^^'AA' = 0 and (p^d^)A^ = 0. (14) 

The solutions to these equations are labeled by additional helicity indices a and a' which 
transform under the little group of the Lorentz group, which is SO{D — 2) in our case. 
Note that in ZA > 4 dimensions helicity of a massless particle is no longer conserved and 
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transforms according to the little group similarly to helicity of a massive particle m. D = A. 
One has 


^A'a' = 0, = 0, (15) 

and for their conjugates 

= 0, (p„5»‘)Af = 0. (16) 

It is always possible to take the solutions to these equations A^/a/(p), A^“(p) (and their 
conjugates) in such a way that 


^A®“(p)Af(p) =p„((T»‘)“, 5^ABv(p)Ayp) =p,.(5'‘)bm- 


(17) 


One can choose the polarization vectors of gluons in the form 

e;,(p|,) . ^ /.(p)(a>BpV(,) ^ ^ V(p)(B-a-)A“'(p) 


(m) 


(pg) 


(18) 


The polarization vectors for massless fermions can also be chosen as spinors while the 
polarization vectors for scalars are trivial. The dependence of the polarization vectors on 
the auxilary light like momenta q reflects the gauge ambiguity in the choice of polarization 
vectors. The dependence on q is always canceled in the hnal gauge invariant amplitude. 

This way one can always write down the scattering amplitude in the gauge theory in 
arbitrary even dimension as a function of the Lorentz invariant products of momenta and 
polarization vectors in terms of the spinor products. 

Concluding this section we would like to comment on the difference between the spinor 
representation of the amplitudes in H = 4, 6 and D = 8,10 dimensions. In H = 4 the 
little group is SO{2) ~ U{1), so its action on the spinors is just a multiplication by a 
complex number. The condition = 0 for p"" = p^(cT^)"" is equivalent to det{p) = 0, so 
the following equality holds: 


= A“A“,. (19) 

On the other hand, one can always use the solutions of the Dirac equation A"(p), A“(p) 
to write 


A"(p)A»(p) = p“ 


( 20 ) 


This means that relation (19) works in ’’both directions”. For given spinors A“ and A“ 
there is always (complex) momentum p such that p“" = A"A" and vice versa one can 
always decompose the light-like momentum p into a pair of spinors using the solution of 
the Dirac equation. This is possible because the product of two D = A Weyl spinors A" 
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and A“ contains 2 x 2 — 1 = 3 (—1 is dne to the little gronp U{1) invariance of A"A") 
independent components, jnst as the light-like momentum p^. 

The same situation occurs in iA = 6. One has 4x2 components in the spinor product, 
and taking into account the action of the SU{2) little group (S'0(6 —2) ~ SU{2) x SU{2)) 
one gets 4 x 2 — 3 = 5 degrees of freedom, exactly as for the D = 6 massless momentum 
P/.- 

For D > 6 this is no longer the case. While it is still possible for a given p^ to find 
solutions of the Dirac equation such that ( |I^ holds, for a given set of spinors in ZA > 6 
there is no unique satisfying this equation [TJ E]. In other words, one may say that for 
D > 6 the spinors obey the nonlinear relations (constraints) [6] and one way to solve these 
constraints is to require that they satisfy the Dirac equation for some light like momenta 
Pf, 0. This means that the spinor helicity formalism in ZA = 8,10 may not be optimal 
(in terms of simplicity) representation or amplitudes. 

3 The on-shell momentum superspace in various 
dimensions 

In the next sections, we discuss the essential details regarding the on shell momentum 
superspace constructions in ZA = 6,8,10 dimensions. 

3.1 D = 6 A/'= (1,1) SYM 

The usage of the on shell momentum superspace allows one to obtain a compact represen¬ 
tation for the amplitudes in supersymmetric gauge theories, which is very convenient in 
the unitarity based computations [T]. 

Consider now the essential part of the ZA = 6 Y = (1,1) on-shell momentum super¬ 
space formalism. The on-shell Y = (1,1) superspace for ZA = 6 SYM was first formulated 
in [laEiin!. It can be parameterized by the following set of coordinates: 

Y = (1,1) D=6 on-shell superspace = {Af, A^, (21) 

where pf and pf are the Grassmannian coordinates, Z = 1,2 and I' = 1', 2' are the 
SU{2)r X SU{2)r R-symmetry indices. Note that this superspace is not chiral. One has 
two types of the supercharges and with the commutation relations 

{?A7'YbJ'} = PAB^rj', 

= 0 . ( 22 ) 

The creation/annihilation operators of the on shell states from the Y = (1,1) super- 
multiplet are 

{Yd, 
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which corresponds to the physical polarizations of the gluon \Aaa)^ two fermions °‘) 

and two complex scalars ) (antisymmetric with respect to I,!')- This multiplet is 
CPT self-conjugated. However, to combine all the on-shell states in one superstate IH) 
by analogy with the TV = 4 D = 4 SYM [20], one has to perform a truncation of the 
full M = (1,1) on-shell superspace [T7] in contrast to the former case. Indeed, if one 
expands any function X (or |r2) superstate) defined on the full on-shell superspace in 
Grassmannian variables, one encounters terms like ~ Since there are no such 

bosonic states A\ ““ in the M = (1,1) SYM supermultiplet, one needs to eliminate these 
terms by imposing constraints on X, i.e., to truncate the full on-shell superspace. If one 
wishes to use the little group indices to label the on-shell states, the truncation has to be 
done with respect to the R symmetry indices. This can be done by consistently using the 
harmonic superspace techniques HZI. 

Defining the harmonic variables uj and which parameterize the double coset 
space 


SU{2 )r SU{2)n 
U{1) U{1) 

we express the projected supercharges, the Grassmannian coordinates 


q 




T] 


T 

a 


Uiq^\ Qa = 


and the creation/annihilation operators of the on-shell states 

0 —, 0 -+, 0 +-, 0 ++, 

^+0 

A ad 


(23) 


(24) 


(25) 


in terms of the new harmonic variables. 

In what follows we consider only the objects that depend on the set of variables 
which parameterize the subspace (’’analytic superspace”) of the full Af = (1,1) on-shell 
superspace 

Af = (1,1) D=6 on-shell harmonic superspace = {A^, A^, fjA }. (26) 

The projected supercharges and momentum generators acting on the analytic superspace 
for the n-particle case can be explicitly written as: 

n n n 

pAB q-^ = ^A^(077-“, q\ = J2Ki'i)vti- 
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(27) 




Now one can combine all the on-shell state creation/annihilation operators (25) into one 
superstate | 12 j) = f 2 i| 0 ) (here i labels the momenta carried by the state): 


la) = {(l>r + (lyriv~v~)i + (l>i~iTT)i + (l>riv~v~UTT)i 
+ {Av~v^)i}\o), 

J_ y-ajay ^ 

^ I aja' 


(28) 


where {XY)i = X- Yi ala- Hereafter we will drop the ± labels for simplicity. As in 
D = A case, we can formally write the colour ordered amplitude as 


A({A^Ai,,„,^i}) = { 0 |nS 2 iS| 0 ), 


(29) 


2 = 1 


Here S is the S-matix operator of the theory, the average (0|... |0) is understood with 
respect to some particular (for example, component) formulation of the theory. The 
invariance with respect to translations and supersymmetry transformations requires that 


= q^An = QaAti = 0 . 


Thus, the superamplitude should have the form: 

Al„({A^,A^.r;„.^i}) = '5‘(p''®)i5'‘(<;'’)'5‘‘(5/i)'P»({A^.A(\.r(„,f;J). 


(30) 


(31) 


where Vn is a polynomial with respect to p and f) of degree 2n — 8 . Note that since there is 
no helicity as a conserved quantum number, there are no closed subsets of MHV, NMHV, 
etc. amplitudes in contrast to the D = A case. 

The Grassmannian delta functions and dehned in this case as 


S^{q^) = ■^_€ABCDSiq^)5{q^)5{q^)5{q^), 

= ^(^^^^^kqA)kqB)kqc)HqD)- 


(32) 


Here the delta function S{Xy is the usual Grassmannian delta function dehned as 6^{X^) = 


X^, where / is the R-symmetry index. In harmonic formulation we simply have 
^(X) = X. 

To extract the ordinary component amplitudes from this supersymmetric expression, 
one has to apply the projection operators. The projection operators are the derivatives 
with respect to an appropriate number of Grassmannian variables. Their explicit form can 
be read from (39). For instance, the projection operator for the i-th gluon is d/dp~d/drjf. 

Gonsider now the four-point amplitude. The degree of the Grassmannian polynomial 
P 4 is 2 ? 7 , — 8 = 0, so Va is a function of bosonic variables {A^, A^} only, just as in the 
D = A case 


M{KA%Va.Va}) = 5\p^^)5\q^)nqA)VA{{xt>^\}) 


(33) 
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At the tree level Vi can be found from the explicit expression for the 4 gluon amplitude 
in components [HI [1^ obtained with the help of the six dimensional version of the BCFW 
recurrence relation mm- Comparing the component expression with and expanding 
(33) in powers of rj one concludes that: ~ 1/st, where s and t are the standard 

Mandelstam variables. So at the tree level the 4-point superamplitude can be written as: 


Af = (34) 

st 

Note that already at the tree level the 5-point amplitude is not so simple mm- Similarly 
to the D = A case, one can obtain the expression up to three loops using the iterated two 
particle cuts. To perform this computation, the following formula for the Grassmannian 
integration is useful: (f f 


+ Af - Qb) 

= {2\)H{h, hf5\qt + (35) 


3.2 L) = 8 A/" = 1 SYM 

Consider now the D = 8 A/" = 1 case. The details of the on-shell M = 1 superspace 
for D = 8 SYM can be found in [15]. It can be parameterized by the following set of 
coordinates 


M = 1 D =8 on-shell superspace = {A^“, A^ , //a}, (36) 

where r]°‘ are the Grassmannian coordinates, a is the little group SO{6) index, and A and 
A' are the spin{SO{7,l)) indices. The R-symmetry group here is U{1)b and 77 “ carries 
the -|-1 charge with respect to U{1)b. Note that this superspace is chiral. 

The commutation relations for the supercharges have the usual form 

(37) 


where the supercharges in the on-shell momentum superspace representation for the n- 
particle case are 


= 


' = 5^A^“(*)Af (i), ^ A"^“(i)r7„, g^' = ^Af (i) 


A. 

dpa 


(38) 


Z =1 2=1 2=1 

The creation/annihilation operator states in the Af = 1 D = 8 on-shell supermultiplet are 


4 >, 0 }, 
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which corresponds to the physical polarizations of the gluon two fermions 

I'ha) and two scalars |0), |0). One can combine them into one ’’superstate” lOj) similar 
to the D = A case 


la) = (a + + l^VaVbAf + |0). (39) 

Here is the absolutely antisymmetric tensor associated with the little group 5'0(6) = 
SU{A). Using the arguments identical to the D = 4,6 cases we conclude that the colour 
ordered superamplitude should have the form: 

A({A-‘‘“,Af= A»(p"‘«')A®(/)P„({A-<“,Af,rfc}), (40) 

where Vn is a polynomial with respect to rj and f] of degree 2n — 8 . The Grassmannian 
delta function is dehned in this case as: 

2=1 

Here is the absolutely antisymmetric tensor associated with the spin{SO{7, 1)). 

Consider the four-point amplitude. The degree of Grassmannian polynomial V 4 is 
2 n — 8 = 0, so as in the previous cases V 4 is a function of bosonic variables and one can 
again write the four-point amplitude in the form 

A({A^“, Af = S\p^^')S\q^)V 4 {{X^^,Xf}). (42) 

At the tree level V 4 can be found from a comparison with the explicit expression in the 
components obtained as the held theory limit of the string scattering amplitude. As in 
the D = 6 case, one has ~ 1/st. So at the tree level the 4-point superamplitude can 
again be written as: 

Af^ = 5%p^^)5^{q^)\. (43) 

st 

Similarly to the D = 4,6 cases, one can obtain the expression up to three loops using the 
iterated two particle cuts. To perform this computation, the following formula for the 
Grassmannian integration is useful, which is similar to the D = 4 case: 

J d\,d\,6^{Xf^'^Pa,h + Xf^Pa^h + <it)^^{Xf^da,h + Xf^^Pa,h “ 9 ^) 

= {A)H{hhf5\qt + q^). (44) 
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3.3 D = 10 M=l SYM 


The D = 10 A/” = 1 SYM supermultiplet of on-shell states consists of the physical polar¬ 
izations of the gluon and the fermion helds. In this case the on-shell momentum 
superspace formalism is not known. The problem is that there are too many r] variables 
[6] (we need 4 rj variables to accommodate all 2^ states in the theory, but the smallest 
representation of the little group SO (8) gives 8). 

However, one can use the indirect symmetry arguments (see the next section) to show 
that the ratio of in H = 10 Y = 1 SYM has the form similar to that in the 

H = 4, 6, 8 SYM theories. One can also use an alternative formulation of the amplitudes 
in the D = 10 Af = 1 SYM theory based on the pure spinor formalism [18] to show 
that at one and two loops the integrand of the ratio of in H = 10 A/" = 1 

SYM indeed has the form identical to that in the D = A case (also component unitarity 
based computations are available up to hve loops [IS]). This strongly supports the above 
mentioned claim. 


4 The A 4 amplitude in SYM theories in various di¬ 
mensions 

4.1 Prom physical to colour ordered partial amplitudes 

The aim is to calculate the multiparticle amplitudes on mass shell. For this purpose, we 
hrst perform the color decomposition extracting the color ordered partial amplitude [T]. 
The relations between physical and colour ordered amplitudes look like: 

.. .p^) = ^ Tr[cr(T“i ... T“’^)]M„(a(p^ .. .p^)) + 0(1/Y,). (45) 

cr£SnlZ„ 

The colour ordered amplitude An is evaluated in the planar limit which corresponds to 
Nc -)■ 00 , -)■ 0 and g^M^c - fixed. 

For the four-point amplitudes the colour decomposition reduces to 

2,3,4)= TUf ^(1, 2,3,4)+ T^Af \l, 2,4, 3) + T^Af\l, 4,2,3) (46) 

where T* denote the trace combinations of SU{Nc) generators in the fundamental repre¬ 
sentation 


rjn2 

^3 


^J,^^a\rpa2rpa4,rpa-i^ _j_ rp^^rpairpa3rpa4,rpa2^ 


(47) 
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What is essential and becomes obvious using the superspace formalism, the four point 
tree-level amplitude always factorizes so that the colour decomposed L-loop amplitude 
can be written in the form 

Af\l, 2 , 3 , 4) = 4°) ( 1 , 2 , 3 , 4)m4(1, 2 , 3 , 4) = 2 , 3 , 4)Mf ^([1 + 2 ]^, [2 + 3 ]") 


or using the standard Mandelstam variables 

4^4, 2, 3,4) = 4°4,2, 3, 4)m4 (s, t) 

It is this M^\s,t) factorized amplitude that we calculate in this paper. 
In general has the form 

M^\sA) = (—coefj X Masterintegralj, 


(48) 


(49) 


where coeh are some monomials of s and f, the Masterintegralj is one of 

the master integrals in the D-dimensional Minkowski space to be evaluated. 


4.2 Dual conformal invariance and the universal expansion 

The D = A J\f = A SYM planar S-matrix in addition to the PSU{2, 2|4) (super)conformal 
symmetry has a new type of symmetry, namely the dual (super)conformal symmetry. One 
can think of this symmetry as of (super)conformal transformations acting on the new dual 
variables and their fermionic counterparts. The dual variables are defined in 
D = 4, 6,8,10 dimensions as 


pr 


= X, 


„AB' 


X. 


^AB' 
i+1 ■ 


(50) 


The explicit form of the generators of the dual (super) conformal transformations in D = 4 
as well as details specihc for D = A Af = A SYM can be found in [20] . 

The dual (super) conformal symmetry is exact at the tree level for a general kinematical 
conhguration. For the amplitudes at the loop level the dual conformal symmetry is, in 
general, broken due to the presence of the IR divergences (see the details, for example, 
in [3]). However, this symmetry is still exact if one considers not the loop amplitudes 
themselves but rather their integrands. For the four-point amplitude this statement mani¬ 
fests itself in the fact that the ratio 4^V4°^ is given by the linear combination of 
the so-called dual (pseudo)conformal integrals with the coefficients given by the rational 
functions of the Mandelstam variables. 

Remarkable that one can dehne, at least, the bosonic part of the dual conformal 
transformations for the D = 10 A/" = 1 SYM [6]. Moreover, it is claimed that the tree 
level S-matrix of the D = 10 A/" = 1 SYM is covariant with respect to the dual conformal 
transformations similar to the D = A case. This would immediately imply the dual 
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conformal covariance of the tree level S-matrix for the D = 4, 6, 8 SYM theories since 
they can be obtained by dimensional reduction |6l UHl |2T] from the D = 10 A/" = 1 SYM. 
This would also imply that the structure of unitarity cuts in all SYM theories is identical. 
Combining all these statements together one can conclude that the integrands of the ratio 
in the D = 6 , 8 ,10 SYM theories have the form identical to those of the D = A 
SYM theory. 

Separate investigation of the dual conformal invariance in the D = 6 M = (1,1) SYM 
fn\ as well as explicit computations of the four-point amplitudes up to 3 loops and the 
spinor formalism based results up to 2 loops in H = 10 [IS] (see also discussion of the 
Y = 1 D = 10 unitarity cuts for the 5 loop integrals in [19]) strongly support that the 
integrands in the D = 4,6,8,10 dimensional SYM theories are indeed identical. 

Consider as an example the four point amplitude m D = A, 6, 8,10 SYM theory. At 
the one loop the contribution to A^l'*/A^^'^ is given by the box diagram with the integrand 


Box Integrand 


st 

F(/c - piYik P2Y{k + P2+ Psy' 


(51) 


Introducing the dual coordinates Xi as (here we omit the spinor indices for simplicity) 


Pi = Xi2, P2 = X 23 , P 3 = X34, Pa = X41, k = Xi3, 
we can rewrite the integrand of this box diagram as follows: 

Boxintegrand = ^ ^ r, . (52) 

^ ' 

X 15 X 25 X 35 XA 5 

This expression is invariant under the dual conformal transformations, except inversions 
and transforms covariantly under the dual conformal inversions I with the weight {xfy 
which is absorbed into the integration measure in H = 4. In higher dimensions this factor 
is left behind, however, it does not depend on external coordinates. The requirement of the 
dual conformal invariance is enough to uniquely fix the combination x\.^x\^/x\^X 2 ^x‘l^x\^, 
i.e. the integrand written in terms of the Box integral 1 /x\^x\^x\^x\^ with the hxed 
coefficient x\ 2 ,x\ 4 - Similar logic is also true in higher orders of PT. 

In combination with the unitarity based method the dual conformal invariance allows 
one to explicitly express the amplitude terms of some dual (pseudo)conformal 

integrals with the known coefficients up to 6 loops. The all loop expression known as the 
BDS ansatz for Aa /was also obtained [22] . 

We base our calculation of the four-point amplitude in various dimensions on a uni¬ 
versal expansion which, as it was mentioned above, is valid for any D. The difference is 
only the dimension of the integration while the integrands stay universal. 
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The expansion for the ratio ^ 4 /^® up to 3 loops is schematically presented below 


= 1 -g St 


+ 9“^ s^t 


-g® sh 


*A\.a 


= i + EMf)(s,t) = 


+ st^ 


+ 2s^t 


+ 2 st^ 


+ 



1 


K 

-f- 


+ st^ 


(53) 


where the connected strokes on the lines mean the square of the flowing momentum. In 
what follows, we consider the D=6, 8 and 10 cases. 


5 Calculation of Integrals 

5.1 Analytical evaluation 


Due to eq.(53) the problem of calculation of divergences is reduced to the scalar master 


integrals, which are universal for any dimension. To evaluate them, we use the dimensional 
regularization. Throughout the paper we accept the following dehnition of the L-loop 
master integrals: 


Masterintegralj = 


ITT 


D/2 


d^ki...d^kL 


D, 


Nurrii. 

Derii- 


(54) 


Since we are interested only in the leading divergences (the leading poles in dimensional 
regularization), the task is essentially simplihed. One has to admit that to calculate the 
leading pole, there is no need to calculate the multiloop diagram itself. The leading pole 
follows from the lowest order singularity due to the nature of the 7^-operation. It is 
valid in any local held theory and guarantees the locality of divergences if the lower order 
counterterms are taken into account. 

Let us briehy recall the main notions of the 7?.-operation |23l [21]. Being applied to 
any Green function T (or any particular graph G, as in our case) it subtracts all the 
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UV divergences including those of divergent subgraphs and leaves the finite expression. 
The use of the 7^-operation is equivalent to addition of the counter terms to the initial 
Lagrangian. The TZ operation can be written in terms of the subtraction operators in the 
factorized form 

no = JJ(1 - K^)G, (55) 

where the subtraction operator subtracts the UV divergence of a given subgraph 7 
(for the minimal subtraction scheme the operator K singles out the l/e"" terms) and the 
product goes over all divergent subgraphs including the graph itself. 

It is useful to dehne also the incomplete 71 operation denoted by 71' which subtracts 
only the subdivergences of the graph G. The full 7Z operation is then dehned as 

7^G = (1 - K)'R!G, (56) 

where K without subscript is Kg- The KTZ'G is the counter term corresponding to 
the graph G. Each counter term contains only the superhcial divergence and is local in 
coordinate space (in our case it must be a polynomial of external momenta). 

The 7Z' operation for any graph G can be defined by the forest formula, but for our 
calculations it is more convenient to use the recursive definition via the 7Z' operation for 
divergent subgraphs (for details and examples see chapter 3 in [25]): 


= 1 - {Kn'y) - ... j G, (57) 

The sum goes over all IPI, UV-divergent subgraphs of the given diagram and the multiple 
sums include only the non-intersecting subgraphs. And 


{KTZ')y G = KTZ'{-f) * G/7, 


(58) 


where G /7 denotes the graph G with the subgraph 7 shrinked to point, the * operation 
inserts the subgraph’s counter term (in our case it is a polynomial of momenta external 
to a given subgraph) into the remaining graph G/ 7 . 

When applying this formula to the diagrams at hand one hnds out that for the n-loop 
diagram the 77'-operation results in the series of terms (we consider only the leading pole) 




+ ... + 


Ai(/i' 


2\e 


(59) 


where the term like ^ comes from the /c-loop graph which survives after subtraction 
of the (n — /c)-loop counterterm. The full expression (59) has to be local, i.e. should 
not contain terms like (log/i)^/e™' for all k,m > 0 while being expanded over e. (For 
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simplicity hereafter we put = /r.) This requirement gives us n — 1 equations for the 
coefficients Ai. Solving them one gets 

-4„ = (-ir+‘—. (60) 

n 

In the case when the hrst divergence appears at k loops (as in the D=6 case where k=3) 
this formula is slightly modihed and looks like 


A, = - ’-A -^ 

n\n — 2)! 


(61) 


This means that performing the T^'-operation one can take care only of the one loop 
diagrams surviving after contraction and get the desired pole term via eq.(61). This 
observation drastically simplihes the calculation of the leading pole. 

To demonstrate how this technique works, we calculate the leading pole of the 5-loop 
diagram in D=6 (see Appendix A). For pedagogical purposes we describe hrst the 
calculation using the full 7?.'-operation and then show how the truncated version using 
eq.(61) works. At the beginning, we dehne the inner momenta, as shown in Figj^ The 
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Figure 1: The choice of the inner momenta 


slashes on the lines correspond to the numerator 

Num = (p 2 -l + P 3 )^(pi + 


(62) 


It is useful to rewrite the hrst bracket as (p 2 — l + Ps)^ = {P 2 — lY + ‘^Pzij >2 — 0 +P 3 - Then, 
having in mind that = 0 we have two terms. 

In the hrst term the numerator {p 2 — lY cancels one of the propagators. The TZ'- 
operation for this diagram is shown in Figj^ There are two IPI divergent subgraphs, the 
three-loop and the four-loop ones. Since one subgraph is inside the other the 7?.'-operation 


contains only the hrst sum in eq.(57). For the graphs framed in dashed boxes one has to 
take the ICTZ' expression. 

We start from the 3-loop subgraph. Since the singular part of the tennis-court 3- 
loop graph equals — l/6e and does not depend on momenta, one actually has to calculate 
the remaining 2-loop graph. To do this, we notice that the one-loop bubble in D=6 is 
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Figure 2: The 7?.'-operation for the hrst term 


proportional to the ingoing momentum squared, so when substituting it into the 2-loop 
graph one propagator will be canceled and the resulting graph will take the form of a 
bubble with ingoing momentum equal to p 2 - Since this bubble is also proportional to the 
square of the ingoing momentum and = 0, this leads the contribution equal to zero. 

We now repeat the same procedure for the 4-loop subgraph. To calculate the ICTZ'- 
operation for it, we have again to substitute the 3-loop tennis-court graph, shrink it to a 
point to get the bubble, which again is proportional to the incoming momentum squared. 
This momentum is not on shell but being substituted into the remaining triangle in Figj^ 
cancels one propagator. This converts the triangle into the bubble which for the same 
reason as above is equal to zero. 

Thus, our conclusion is that the hrst term gives zero contribution. 

Consider now the second term. The T^'-operation in this case is shown in Figj^ Again 
we have two divergent subgraphs, one inside the other, which give two contributions. The 
hrst one contains the 2-loop graph. Its singular part can be calculated via the T^'-operation 
shown in Figj^ The divergent subgraph here is a triangle 

/ = -^2!>3(2p, -k+p,) + finue term. (63) 


Substituting it into the remaining triangle obtained by shrinking the hrst one to a point 
we get the integral 


r 2p^{2p2 - k + pi) (f 

J {p2-ky{ky{pi +ky 


1 

6e 


2p3(5p2 + 4pi) -|- finite terms 


Now we have for the 2-loop graph 

2p3(5p2+ 4pi)/i" 
■ e2 36e2 


(64) 


( 65 ) 
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Figure 3: The T^'-operation for the second term. The dot corresponds to the numerator 
2P3(P2 - 1 ) 


R' 



Figure 4: The 7?.'-operation for the 2-loop subgraph 




From this equation requiring the cancellation of the log/r/e term we get 

^ 5t + Au t — As 
" 72e2 72e2 

Thus, the hrst contribution to the 77'-operation is 

/ 1 \ (f — {t — 4s)/i^'^ 

V ^) 72^2 - 216-263 


( 66 ) 


(67) 


For the second contribution the problem is reduced to the 4-loop counterterm. To 
compute it we again use the the 77.'-operation. The corresponding graphs are shown in 
Figj^ The calculation repeats the one performed above (63) and one gets the 77'-operation 
for the 4-loop subgraph 


4e 


n' : A.At 


■i) (-^2 p3(2p.-. + p.) 


This gives for 


Aa — 


2p3(2p2 -k + pi) 
4-36 


( 68 ) 


(69) 
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Figure 5: The T^'-operation for the 4-loop subgraph 


Thus, the JCTZ' for the 4-loop subgraph is 

^ 2p3(2p2 ~ ^ + _ 2p3{2p2 - k + Pi) ^ g 2p3(2p2 -k + pi) 


4 ■ 36e2 


36e2 


4 ■ 36e2 


Next, we substitute the obtained coefficient into the second part of the T^'-operation for 
the 5-loop subgraph Figj^ The resulting triangle integral including the numerator has 
the form (64) which gives 


3/i^ 


4 ■ 36e2 




(71) 


Summing up eqs.( 67||7l ) for the full 7?.'-operation for the 5-loop graph one has 

^/. -4^ /i2^(t - 4s) /i"3(t-4s) 


2 ■ 216e3 


4■216e3 
^2 1„„2 


(72) 


expanding over e and collecting the terms of log/r/e^ and log /i/e we get two equations 
to determine A^: 


log/i : 57 I 5 + 2/2/216(t - 4s) - 3{t - 4s)/4/216 = 0 (73) 

log^ /i : 257 I 5 + 4/2/216(t - 4s) - 3(t - 4s)/4/216 = 0 (74) 


Solution to these eqs is 


(t — 4s) s —t/4 
20 ■ 216 “ 30-36 


(75) 


Consistency of the two equations serves as a check of correctness of the calculations. 

One can arrive at the same result in a shorter way using the truncated 77'-operation 
and evaluating only the term with the one-loop subgraph (71). Indeed, using relation (61) 
for u = 5, fc = 3 one gets 


1 ^ 1 3(t — 4s) s — t/A 

15"^^ “ ~15 4-216 “ 30 - 36 ■ 


This gives us the result listed in Table 


( 76 ) 
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Figure 6: The choice of the inner momenta 


As the second example we consider the evaluation of the the diagram in D=8 by 
means of the truncated 7?.'-operation. Again, we define the inner momenta, as shown in 
Figj^ Keeping in mind that p? = 0 we rewrite the numerator in the following way: 

Num = (pi + nY = 2(pin) + v?. (77) 

The term proportional to cancels the corresponding propagator. Thereby, after shrink¬ 
ing the 2-loop subgraphs we get the bubble diagram which is proportional to the on-shell 
momentum squared and is equal to zero. Thus, the numerator finally gets the form 
Num = 2(pin). 

Consider now the 77' - operation for this diagram having in mind that we are interested 
only in the one loop remaining graphs shrinking the rest of the diagram to a point. Since 
the upper right and left boxes transform into bubbles after shrinking the rest part of the 
diagram and these bubbles are again proportional to the on-shell momentum squared, the 
answer for both of them is zero. Therefore, the only one loop graph that survives is the 
lower box. The 77' - operation is shown in Figj^ 



Figure 7: The 77'-operation for the 3-loop graph. The dot corresponds to the numerator 

2(pin) 


The calculation of the double box subgraph is performed again using the 77'-operation 
shown in Figj^The integral for the left box has the following form: 


/ 


k‘^{n 


kY{l-kY{pi + kY 


-h finite terms. 

6e 


(78) 
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Figure 8 : The T^'-operation for the double box 


Since its singular part does not depend on momenta, the right box is the same. For the 
hrst and second triangles in Figj^the integrals are 

f 2x1^ + 2{pin) 


J k‘^{n — kY{pi + ky 24e 

+ 2 ■ 2{pin) + 2 (p 4 n) + t 


finite terms, 

I / + finite terms. (79) 

J {pi + iy{n - iy{p2 + P3 - ly 24e 

Substituting expressions ( 78|[79 ) into the 77'-operation we have for the 2-loops subgraph 

2n‘^ + 2{pin)\ 1 

-25—jfe 


4 /^^' ef + 2 ■ 2{pin) + 2{p^n) + t\ 1 

I-- 1^-/^ 


Evaluating the coefficient A 2 from the requirement of cancelation of the logp/e term one 
gets 

+ 3 ■ 2{pin) + 2{p4^n) + t 

-• (oJ-J 


A2 — 

This gives the ICTZ' for the double box 
/C77'Double Box = A 2 ^— 


2 ■ 6 ■ 24 

An^ -I- 3 ■ 2{pin) + 2(p4n) -|- t\ /i^ 
2 ■ 6 ■ 24 


An^ -I- 3 ■ 2{pin) -|- 2{p/^n) -|- t 


(82) 


2 ■ 6 ■ 24 e2 

We now turn back to the 77'-operation (Figj^. Substituting (82) and having in mind the 
numerator 2(pin) we hnally get the triangle with 3 powers of internal momentum in the 
numerator. Performing this integration one has 

Ai t {Sy — 2st + s^) 


3!4!5!3e3 


According to eq. (60), this gives 


1 t (Sy — 2st -|- s^) 

^3 = -Ai = — -——- - 

^ 3 3!4!5!9 


(83) 


(84) 


(3) 

This way we get the expression for the leading pole of the diagram . 

Applying the described truncated 77'-operation we calculated all the leading poles for 
the diagrams in D = 6 up to 5-loops and in D = 8,10 up to 4-loops. The results are 
presented in Tables and 
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Table 1: The leading poles of the diagrams up to 4-loops for D = 6,8 and 10 
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Table 2: The leading poles of the diagrams in 5-loops for D = 6 
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5.2 Numerical evaluation 

Since in higher loops the evaluation even of the leading pole happens to be a compli¬ 
cated task, we performed a numerical check of our calculations. To evaluate the graph 
contributions to the leading pole numerically, we hrst go to the Euclidean space and then 
use the alpha-representation and the method of sector decomposition m- 

For the diagrams under consideration the leading poles {L.P.) are polynomials of s 
and t. For the diagram G with the degree of divergence equal to 2N the leading pole is 

N 

L.P.{G{s,t)) = Y,C^,N-^ s* (85) 

i =0 

The coefficients Gij can be calculated performing the differentiation of the integrand over 
s and t ^ 

L.PiG{s,t)) = ^ s* t^-^L.Pi G,,M-^ ), = f fr G{s,t) (86) 

i=0 

For i + k = N the integral for Gi^k becomes logarithmically divergent and L.P.{Gi^k) is a 
constant which we calculate using the sector decomposition technique. It turns out that 
it is more convenient to calculate the graph Gi^k with massive lines and to put s = t = 0 . 
The leading pole of such a graph is exactly the same as for the massless graph 

L.P (G,,k\s,t=o) =L.P{G,,k); (87) 

however, choosing all momenta to be zero one gets additional simplihcations in the sector 
decomposition technique. 

The standard form of the alpha-representation in the Euclidean space for a graph 
without the numerators is 


, TT ^ n 

U)DL/2 roo _y/U-J2m,a, 

I{s,t,mi) = — - / dai...dan e 

n r(A.) to 


( 88 ) 


where A* are the powers of the propagators, the functions U and V are the polynomials 
of the alpha-parameters (a*) of the order L and L -|- 1, respectively. The polynomial V 
linearly depends on the squared combinations of the external momenta of the graph. In 
our case 

V = s ■ Ps{ai) + t ■ Ptitti) . (89) 

In the general case (for the graphs with the numerators), the alpha-representation 
can be obtained from the generating function (for details see [2l]), which allows one to 
construct the alpha-representation for any particular numerators. This procedure though 
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straightforward is rather lengthy. Instead, since we have very specihc numerators, we 
construct the generating function in terms of the so-called dual graphs. (Some examples 
are given in Figsj^and 10). Note that this is always possible for the planar diagrams. In 
the dual representation the numerators can be expressed as additional lines with negative 
powers of the propagators (see the dotted line in FigJIo|). 

The propagators with negative powers can be treated as the ordinary propagators but 
instead of integration over the alpha-parameter one should differentiate with respect to it 
and after that set it to zero 




(vr)W2 

n 

n r(A,) 


n+k 


n (-3».) 


Ki 


y^i=n+l 


n 

f°° dai...dan -VlU-^mja^ 

Jo 


^ n +1 —• • • — —0 


(90) 

Here the parameters a„+i,..., correspond to the numerators, and Kn+i, ■■■, Hn+k are 
the powers of the corresponding numerators. The advantage of this approach with respect 
to the standard one is that we have only two polynomials U and V and no additional 
terms appear. 

For the graphs without the numerators the monomials in U are generated by the 1- 
trees and have a degree of L. The momenta corresponding to the alpha-parameters in each 
monomial should be linearly independent (assuming that all external momenta and their 
combinations are zero). In terms of the dual graphs they should not form loops. Each 
monomial in V is generated by the 2-trees and contain the factor equal to the squared 
momenta flowing from the one part of the 2-tree to the other. 

Note that for the polynomial U in the dual representation one has to treat all exter¬ 
nal vertices as identical (equivalent to setting all combinations of external momenta to 
zero) and avoid the combinations with chains that contain two external vertices. The 
polynomial V is constructed in the same way, but monomials must contain exactly one 
chain between external vertices, and is multiplied by the combination of external mo¬ 
menta squared flowing through the lines that belong to this chain. Actually, this chain 
works like a cut in the standard (not dual) representation and the momentum factor is 
the momentum flowing through this cut (this is exactly what 2-trees do). Since in our 
case we have only momenta s and t in V, this means that for Pg we need the chains 
that connect the f-vertices and for Pt the chains that connect the s-vertices (see example 
below). 

To illustrate the equivalence of the dual graphs approach with the usual 1- and 2- 
trees, we consider Erst a simple example without the numerators (see Figj^. 

To construct the polynomial U for this graph, we need to find all possible pairs of 
alpha parameters which do not form the loops in the dual graph while all external vertices 
(marked by s and t) are considered as identical. One has the following combinations of 
two and three alpha-parameters which form the loops in the dual graph: 


{q:iQ:2, Q:iQ:6) 0:30:4, 0:30:5, 0:30:5, 0:40:5}. 


(91) 
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Figure 9 : The double box diagram and its dual graph 


{010307, 01O4O7, 01O5O7, 02O3O7, 02O4O7, 04O6O7, 05O6O7, 02O5O7, 03O6O7}. ( 92 ) 

Since the polynomial U has degree of 2 we take all possible pairs of alpha-parameters 
except for those listed in ( 91 ) and get: 


U — O1O3 O1O4 -|- O1O5 -|- O1O7 -|- O2O3 -|- O2O4 -|- O2O5 -\r O2O7 -\- 030g O3O7-I- 


( 93 ) 


-I-O4O6 "F O4O7 -f- OgOg -\- O5O7 -f- 0g07 

It is easy to check that this polynomial is exactly the one constructed using the 1 -trees. 
The excluded combinations of alpha-parameters ( 91 ) are the ones for which the rest of 
the graph does not form the 1-tree. 

To get the polynomial V ( 89 ), one needs to find Pg and Pf. Pg consists of the combina¬ 
tions of three alpha-parameters which contain the chain of lines connecting the t-vertices 
and does not form any other loops in the dual graph (the same for Pt but the chain should 
connect the s-vertices). There are only 4 chains that connect the t-vertices: 

{0(20:6, 0(30(5, 0(20(50(7, 0(30(50(7}. ( 94 ) 

Each monomial in Pg is the combination of three alpha-parameters which contain the 
chains from ( 94 ) and do not form any additional loops ( 91 ),( 92 ). Thus, for Pg we have: 

Pg = 0(20(6(0(3 -F 0(4 -f- 0(5 -|- 0(7) -|- 0305(0(4 -f- 0(2 + 0(6 -|- 0(7) -|- O2O5O7 -|- O3O6O7. ( 95 ) 

For Pt there is only one chain ({0(40(40(7}) that connects the s-vertices, so Pt has form: 

Pt = O1O4O7 ( 96 ) 

This concludes the construction of all ingredients needed for the alpha-representation of 
the double box graph. 
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Figure 10 : The tennis-court diagram and its dual graph 


Consider now the ’’tennis-court” diagram (Fig. 10 ). The dashed line (an) in the hgure 


corresponds to the numerator and in terms of the dual graph is treated as a normal line 
with the propagator in the negative power. The polynomial U contains the monomials 
with all possible combinations of 3 alpha-parameters except for the ones which form the 
loops in the dual graph. For this graph, the loops are formed by the following combinations 
of alpha-parameters: aia2, a3a4, asag, aiay, a2a7, aian, a2aii, ayan and a2a3a8, 
a4a8aii, a4a5aii, etc. (totally there are 21 combinations which contain 3 parameters and 
22 combinations with 4 parameters). This results in 80 terms for U. The polynomial Pg 
contains the chains that connect the t-vertices and contain no other loops, i.e., a2a3a9a7, 
etc., (totally 32 terms). For Pt there are 25 terms. 

Thus, this diagram in the alpha-representation looks like 




dai-.-daio -{sPs+tPt)/u-J2^j(^j 

- p 

JJdl2 


aii=0 


( 97 ) 


In the case of H = 6 , this diagram diverges logarithmically, and we can set both s and t 
to zero and calculate the integral 


Gift ^\s = 0,t = 0,mi) = j 


f°° dai...daio “ F) 


q:ii=0 


( 98 ) 


In the D = 8 case, the graph has degree of divergence equal to 6 . Hence, one has four 
contributions proportional to s^, sP and P. Taking the corresponding derivatives one 
has 


G'tth = 0 ,i = O.m.) = 


oii=0 


( 99 ) 
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Here the extra multiplier {—Psy{—Pt)^~^/U^ originates from differentiation of (97) with 
respect to s and t. 

For the integrals like (98),(99) the leading pole can be extracted using the sector 
decomposition method. To do this, we adopt the Speer-like strategy [27] for the dual 
graphs. For the tennis court diagram this strategy produces 390 sectors, for D = 6 all 
the sectors contribute to the leading pole, but for 77 = 8 the leading pole is present only 
in 144 sectors out of 390. 

We performed the numerical evaluation of the leading pole for the diagrams up to 4 
loops in 77 = 6 and 77 = 8, and up to 3 loops in 77 = 10. The obtained results are 
in good agreement with the analytical values, as can be seen from Table The overall 
computational time is about 300 hours at the 140 core cluster. 


graph 

term 

numerical 

exact 

JW) 

-'i 


0 

0 

j{4) 

-'2 
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43 

7° 

0.0208328(7) 
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-210 



-6.6661(10) 

-20/3 

j(4) 

42 

sH 

0.888900(24) 

8/9 


sH‘^ 

-0.1111105(7) 

-1/9 



-20.4765(8) 

-430/21 

j(4) 

43 

sH 

0.444420(25) 

4/9 


sH^ 

-0.0555541(10) 

-1/18 


Table 3: The numerical values for some sample 4-loop graphs in 77 = 6 (left) and 77 = 8 
(right). The values shown in the right table are multiplied by 6!5!4!3!. 



















6 Summary of the leading pole evaluation in various 
dimensions 

We summarize here the results of calculation of the leading poles in various dimensions. 


6.1 D = 6 A/'= (1,1) SYM 

Summarizing one has for the leading poles (L.P.) [13] 


L.P. = 2stg^ 


9' 


s + t + + + IsH + + t^ 

+ 9 -- +9 - 


6 e 


36e2 


216e3 


( 100 ) 


The leading powers of s ant t remind the geometrical progression while for the mixed 
ones there are too few terms to make any guess. If taking the geometrical progression 
seriously, one gets 


E 

n=l 


gh 

6 e 


6 e 


1 _ 21 ^ 
6 e 


( 101 ) 


which looks precisely like the D=4 Yang-Mills theory with the replacement g^ —)■ g^s and 
in the limit when e —)■ -|-0 tends to —1 when s < 0 and to cx) when s > 0. A natural 


question arises whether one can prove eq.(lOl) to be correct. For this purpose consider 


the sequence of diagrams appearing in the loop expansion (53). 


We start with the inhnite sequence of diagrams originating from the graph ±2 


(3) 


by 

adding the boxes to the left and to the right. This gives us the diagrams 
etc. Performing the 7?.'-operation and looking for the surviving one loop diagrams one can 
notice that they can stand either on the left or right edge of the diagram or in the middle. 
But the tennis court subgraph, l'^\ can be present only once since it is a three-loop block 
and if it stands twice the order of the pole drops by two. And the diagrams containing 
only boxes do not diverge. Hence one is left with one option: the one loop graph stands 
at the edge. This is shown in FigJTTj 


R': 



Figure 11: T^'-operation for the ladder-court graph. Shown are the one loop surviving 
graphs only. The dotted line denotes the contracted graph for which the KTZ' has to be 
taken. 


Consider hrst the case when the tennis court is situated at the edge and the boxes are 
added to one side. Then, since the triangle graph itself is equal to l/2e, the KlZ'Gn-i = 
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(—l)”y4„_i and taking into account eq.(60) one gets 


nAn - 

where n is the total number of loops. It has a solution 

” 2" nV 


(102) 


(103) 


where the constant c can be found from the n = 3 case and is equal to 8. 

For the general case, when the boxes are added from both sides one has to take the 
sum of all diagrams with k boxes to the left and n-2-k boxes to the right and sum over k 
from 0 to n-2. The sum obeys the equation 




(104) 


with the solution 


= (-ir-, 

nl 


c = 2. 


r(4) r(4) A5) A5) 

2 l-'g ,±2 


(105) 

are reproduced by 


One can check that the corresponding diagrams 
eqs.( fl03p05| . 

Note that factorial comes inevitably due to the linear nature of eqs.( 102[l04 ). It 
reflects the fact that the boxes themselves are hnite. Therefore, this type of behavior will 
take place for all series of diagrams in this theory, in particular the ones that start with 
the diagram and the following ones. In any new order of PT the new series starts. 

Knowing the n-th order coefficient one may sum all of them and get 


n=3 


g^s 


n—2 


n\ 


2 e^ 

^4^2 


g^s 


,9^s. 


-A 


oo 


s < 0 
s > 0 


when e —)• +0 
(106) 

One can see that we get the same result in asymptotic as in eq.(lOl), but we have taken 


into account only one sequence of diagrams which give the main contribution to the 
coefficients. 

One can continue this procedure and sum the diagrams of the next series which starts 


(4) 


The difference here is that is not a constant but is 


with the 4-loop diagram 
proportional to t — s. As for s, it stands outside the integration and is not changed but 


t is replaced by t' corresponding to the contracted diagram in analogy with FigTl and 
has to be integrated over triangle giving both s and t. Thus, we have two relations, one 
proportional to s, and the other to t. Considering the sequence of diagrams where the 
boxes are added from one side, one gets the relations 


nAi = -Q^h-i, 




(107) 
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Solution to these relations is 


4* — 


- 1 )’" 1 


^n—3 




1 (- 1 )^ 1 


n\ 


2 6 


n—3 


n\ 


2 2 


n—3 


n\ 


(108) 


One can check that eq.(108) is valid for the diagrams and Summing up the 

diagrams where the boxes are added from both sides, similarly to the previous case, one 
has 

= ~^n-l + (109) 


with the solution 


Ei = 


(- 1 )^ 1 
3’"-3 n\' 


E- ^ 


2 3’^-3 


n\ 


2 ^ ^ n! 


( 110 ) 


These relations reproduce the diagrams and 

Having these coefficients one can again calculate the whole series 


3-E 


2 \ n—2 

9 s\ 1 


3e 


t27e^ 
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i 3 
's 2 I 


OO 


sg^s^ 


\[i + m)] .<0 


s > 0 


, g'^s l.g'^s.n l.g‘^s.n 

£3. - 1 - ^-(—)-( —) 

3e 2^3e^ 6^ 3e ^ 


when e —>■ +0, 


( 111 ) 


n=4 




g'^s 


2 c/ 4 s 2 




2 ^ e ^ 


\(^f 


-)■ 


27 


OO 


g^S _ 1 .£^ x 2 _ 1^£!^')3 
3e 2 ^ 3e ^ 6 ^ 3e ^ 


[i + l(t)] + i[i + i(^)] = -i ^<0 


s > 0 


when e —)■ +0. 


( 112 ) 


Thus, we see that in the limit e —)• +0 when s < 0 the hrst series (106) tends to a 


constant and the second to a constant plus the hrst pole (111 [112 ). Obviously, the third 
series, which starts from 5 loops, will tend to a constant, the hrst pole, the second pole 
and so on. This new series has to be summed again. It has a feature common to all the 
sequences, namely it falls as 1/n!. At the same time the number of diagrams in each order 
is expected to be of the order of n!n^a”. For s > 0 all the series diverge when e —>■ +0. 

The same is true to the diagrams in the t-channel with the obvious replacement s ■H-1. 
In the next section, we will show how these results can be promoted for the general case. 
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6.2 D = 8M=1 SYM 


Summarizing one has for the leading poles 
L.P. = —st 

« 1 16770 s 6 - b?>QsH + 4125 ^^^ 
+ 9^ 


2 ^ , + 

^ 3!e ^ 3!4!e2 


g4 15s^ — S^t + §2^2 _ _|_ 

3 3!4!5!e3 


(113) 


3845^^3 + 4125^^4 - 536st3 + 16770f® 


63 


3!4!5!6!e4 


This expression does not look simple even for the leading powers of s or f though numer¬ 
ically one almost has a geometrical progression with slightly rising coefficients. 

Consider now the inhnite sequence of horizontal boxes and apply the 77.'-operation. 
The difference from the previous case is that now the box diagram is divergent, and 
performing the 77'-operation and looking for the surviving one loop diagrams one has to 
consider the diagram in the middle as well, 
changes. It is shown in FigJT^ 


Therefore, the graphical form in Fig. 11 



Figure 12: 77'-operation for the ladder type graph. Shown are the 1-loop surviving graphs 
only 


Since the triangle diagram in D=8 is equal to — l/4!/e and the bubble one to 2/5!/e, 
this gives us the recurrence relation for the leading pole terms 


n—2 


2 2 

^ ^ ^k^n—l—k') Ti ^ 3 


(114) 


k=l 


with Ai = 1/31. Starting from this value one can calculate any An though the explicit 
solution is not straightforward. However, since we actually need the sum of the coefficients 
we apply the summation multiplying both sides of eq.(114) by z being je 

(115) 


OO r) °° r) °° n—2 


4 !^ / 5 ! 

n=3 n=3 n=3 k=l 
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Denoting now the snm by and performing the interchange of the order 

of snmmation in the nonlinear term we get 


d 2 2 


(116) 


Having in mind that 


S3 — Si + Aiz — A2 z‘^, S2 — Si + Aiz, Ai — —, A2 — 


3!' 


3!4!’ 


one hnally gets the eqnation for S = S] 


1 2 _ 2 
““3! ■ 


Solntion to this eqnation is 

S(^) = 


4tan[2;/(8\/l5)] 

1 — tan[2;/(8A/l5)]-\/5/3’ 

The expansion of tan z contains the Bernnli nnmbers 

,_i22’"(22” - 1)5; 


(117) 


(118) 


tanz = y~~^(— 1)’’ 

n=l 


(2n)! 


2n 2n-l 

>2 • 


Being snbstitnted into eq.(118) it gives 

E(z) = -(z/6 + 127144 + 1272880 + 7^7414720 + ...) 


(119) 


Remind that here z = This series reproduces the diagrams and /7 

which give the main contribution to the amplitude PT series. 

The function S given by eq.(118) has an inhnite sequence of simple poles and thus has 
no limit when e —?• 0. This is similar to the geometrical progression with non-alternating 
series. 

The next sequence of diagrams comes from the tennis-court one supplemented by 
boxes from both ends. Here one again has the nonlinear terms like in eq.(114) and the set 
of three sums proportional to s^, st and like in eq.(109). This looks more complicated 
but is simplihed by the fact that the tennis court diagram appears only ones, so equations 
are in fact linear and generally look like 


^'court — C — S 


court 


+ E 


court ^box 


( 120 ) 


with Tibox given by eq.(118) above. 

Since Tn,ox has a singular behavior when e —)■ 0, so does T^court- Moreover, for any hnite 
number of sequences of this kind one will always have a singular behavior. Only when one 
has an inhnite number of them, which is actually our case, one may avoid this singularity. 
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6.3 D = 10 M=l SYM 


Summarizing one has for the leading poles 


L.P. = —st 




+ 9 


, 8 s^ + 2 sH + 2 sP + 8 P 


5!7!e2 


( 121 ) 


@2(20955^ + llSs^t + 33sH^ - llsH^ - llsH^ + 33sH^ + llSsf® + 2095t^) 
^ ^ 5!7!7!45e3 

8 32(2112188805^° + 753490s°f - 1395096s®f2 + 1125763s^t3 _ giegiOs^f^ 

^ ^ 13!7!7!5!5e4 

+843630s^f5 - 916916s^f° + 112576353^^ - 139509652^8 + 7534905t° + 211218880f^°) 

13!7!7!5!5e4 


One can construct the recurrence relations here as well. For the box type diagrams 
one has a relation similar to the D=8 case but since the one loop box has the numerator 
(5 + t) one has two separate expressions like in eq.(107). At the same time they are 
nonlinear like eq.(114). One has 


nAl^ = 
nK = 


with Al = A\ = 
calculated above. 


n—2 




7!' 


3-7! 


k=l 


—A* - -Y 

3 ■ 5! 7! *^-1 


n—2 


7! 


k=l 


E — 4 * 4 * — 4 * 4 ^ _L 

1 ^^k^n-l-k ^k^n-l-k ^k^n-l-k + 


1/5!. These equations also reproduce all the 


( 122 ) 

(123) 

g^k^n-l-k j 

ladder type diagrams 


7 All loop recursive equations 

It is possible to construct the recursive relations for all the diagrams. They actually follow 
from the analysis of the series and the way how the diagrams are constructed and enter 
the 77'-operation. 

Indeed, one particular way to obtain the dual conformal invariant diagram of the n-th 
order from the (n-l)-order is to follow the so-called ’’rung rule” [2H1 122], which states 
that one has to take the diagram of the (n-l)-th order and insert a line connecting each 
pair of the neighboring lines multiplying by a factor equal to the square of the total 
momentum (see Figj^left). This is true for all the dual conformal invariant diagrams 
with triple vertices, which are proportional to the common factor st, and does not include 
the diagrams with numerators without this factor and the diagrams with quartic vertices 
|19j . However, the latter ones do not contain the leading poles. 
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Figure 13: The ’’rung rule” (left) and its application for the box diagram (right) 


To demonstrate how it works, let us take the simplest box diagram. The application 
of the rung rule reproduces the horizontal and vertical double boxes with appropriate 
coefficients (see Figjl^ right). This way, starting from a single box and applying the rung 
rule one adds boxes to the left and right and to the top and bottom, thus creating new 
s-channel and f-channel diagrams. 

What is important, to get the leading pole, the diagram should contain the maximally 
divergent subgraphs which being shirked to a point leave the one loop divergent graph. 
Integrating over this graph gives the desired leading pole. The topology of all the graphs 
constructed via the rung rule leaves just one possibility: the maximally divergent graph 
should have a box at the left (right) or top (bottom) edge. All the other graphs either do 
not contain the maximally divergent subgraphs or, if they do, their shrinking to a point 
leaves the one loop subgraph, which contains a single leg with light like momenta and 
hence is equal to zero. This property singles out the set of maximally divergent graphs. 

Consider now the T^'-operation for this set of graphs. For the the s—channel type 
diagram it has the form graphically presented in Figsj^and 12 where now the (n-l)-th 
order diagrams in dotted boxes contain both the s-channel and f-channel contributions. 
Here we use the above mentioned property that the diagram of interest always contains 
the box either on the left or on the right edge. The K.TZ'Gn-i is a polynomial in s and f; 
however, s is a common factor while t for Gn-i contains the integration momentum over 
the last loop. Substituting the explicit form of s and t and integrating over the triangle 
by introducing the Feynman parameters one gets the desired recursive relation. 


7.1 D = 6 A/'= (1,1) SYM 

In the case of H = 6 it looks like 


1 pX 

dx dy {Sn-iis,t')+Tn-i{s,t')), n>4 (124) 

Jo 

where t' = tx + uy, u = —t — s, and S 3 = —s/3, T 3 = —f/3. Here we denote by S'„(s,f) 
and Tn{s,t) the sum of all contributions in the n-th order of FT in s and t channels, 
respectively. The same recursive relation is valid for the f-channel diagrams with the 
obvious replacement s -H- f. Due to the s — t symmetry of the amplitude, one should have 


nSn{s, t) = — 2 s 
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Tn{s,t) = Sn{t,s). The coefficient An{s,t) of the n’th order pole is the sum 

A„(s, t) = Sn{s, t) + Tn{s, t). 


Eqs. (|124|) reproduces all the diagrams calculated above and may serve as a generating 
LS in all ord( 

S'4(s,t) = 


function for the diagrams in all orders. The hrst few terms are 

2 s^ + st 


S5{s,t) = 

Se{s,t) = 


36 ’ 

5s^ + 

540 ’ 

25s^ + 5s^t — 


(125) 


19440 

The Mathematica hie with a simple code, which generates the Sn and T„ polynomials up 
to given order is available. 

Since we know now all the leading pole contributions in all orders of PT, it is tempting 
to sum them over. However, explicit solution of the recursive relation (124) is problematic. 
Instead, we proceed in the following way: we multiply both sides of eq.(124) by 
and sum over n from 4 to oo, 2 ; being je. Then on the left hand side one has a derivative 


= 2 s 


dx dy {E 3 {s,t',z)+ j:s{t',s,z))\ 

t'=xt-\-yu‘ 


(126) 


Using the fact that ^4(8, t, z) = ^, 3 ( 3 , t, z) + S'3(s, t)z^ one gets the equation 


^^ 3 ( 5 , t, z) = 3 S 3 z‘^ + 2 s 


dx d|/ (S 3 (s,f', 2 ;) + S 3 (T,s, 2 :))| 

t'=xt-\-yu' (127) 


Since the hrst divergence appears only in the third loop order the function of interest 
is S(s,f, 2 :) = = z~‘^T, 3 {s,tt z). Substituting it into eq.(127) and taking 

into account that 83 ( 3 , 1 ) = —s/3 one gets the hnal equation for the sum of PT series 


—T,(s,t,z) = s -S(s, t, z) + 2 s 

dz z 


dx dy {E{s,t',z)+ E{t',s,z))\t'=xt+yu- (128) 


One has the same equation in the t-channel 


—S(t, s, z) = t -S(t, s, z) + 2 t 

dz z 


dx dy {E{s',t,z)+ E{t,s',z))\s'=xs+yu- (129) 


Summing them up one gets for the total sum E(s, t, z) + E(t, s, z) 
d 2 

— (S(s, t, z) + S(t, s, z)) = {s + t) -[S(s, t, z) + S(f, s, z)] 

CLZ Z 

p 1 px 

+ 2 s dx dy [E{s,t',z)+ E{t/s,z)]\t'=xt+yu 


(130) 


+ 2 t dx dy [Tj{s',t,z)+ J:{t,s',z)]\s'=xs+yu- 
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The behavior of the solution to this equation is dehned by the hxed point, i.e. the 
zero of the right hand side. As z —>■ cxd one can neglect the second term and under the 
assumption that the hxed point is a constant get the following conjecture: 

S(s, f, 2 ;) + S(t, s, ; 2 ) = —1. (131) 


Consider now how this hxed point is approached. The sign of the derivative is propor¬ 
tional to s + t = —u. In the case when u < 0 the derivative is positive above the hxed 
point and negative below it. So if the initial value of S is above the hxed point, it will 
increase and if it is below it, it will decrease. This means that the hxed point is unstable. 
On the contrary, if u > 0 the sign is changed and the hxed point is stable, the solution 
tends to it as z —)■ -|-cx) or e —)■ -|-0. Therefore, the stability properties depend on the 
kinematic region. For u > 0 the hxed point is stable and the theory is hnite in the limit 
6 —y TO. 

The situation will be the same for other partial amplitudes. In the (s,u) channel the 
theory is hnite if f > 0 and in the (t,u) channel it is hnite if s > 0. Unfortunately, all 
three conditions are incompatible since s T t T n = 0 and one can not have all of them to 
be positive. 


7.2 L) = 8,10 A/" = 1 SYM 

Consider now the cases of D=8 and D=lt). Here according to FigJT^ one has an addi¬ 
tional nonlinear term, therefore eq.(124) is modihed. Note that for this last term when 
integrating over the loop on both sides one has functions of s and t. Replacing t by t' one 
should have in mind that on the left t' = {I — pi) and on the right t' = [I TP 4 ), where I is 
the integration momentum. This means that while integration one gets the mixed terms 
like This can be taken into account and one gets an equation which looks the 

same way for D=8 and Zi)=10. For D=8 one has 


nSn{s,t) = - 23 “^ dx dyy{l-x) (S'„_i(s,F) + Ai-i(s, 


'0 

n—22k—2 


T s 


/ I n—2 2k —2 

dx x^{l — xY ^ ^ 

k=l p=0 


=1 pKp + 2)! dt'P 


dP 

— (S,{s,t')+T,{s,t'))x 


dP 


X^{Sn-l-k{,S,t') +Tn-l-k{s,t'))\t'=-sx {tSx{l-x)Y, (132) 

where Si = Ti = ^. Equation(|l32|) reproduces all the above calculated diagrams. 
The terms with the derivatives in the second term do not contribute so far. 
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In the case of D=10 one gets analogously 


nSn{s,t) = -s^ / dx dy y‘^{l - xf {Sn-l{s,t') + T„_i(s, t'))|t'=te+yn 
JO JO 
n—2 3k—2 


+ s' 


n—Z 6 K—Z 

' / dx a;^(l — x)^ ^ ^ 

k=l p=0 


=1 p'^n pKp + 3)! dt'P 


1 dP 

^(s,is,t')+ms,t'))x 


dP 


X—{Sn-l-k{s,t') +Tn-l-k{s,t'))\f=_sx (tsx(l - x))^, (133) 


where S'! = |i, Ti = ■^. Here the terms with the derivatives work starting from 3 loops. 

Equations (132), (133) can be summed the same way as in the D=6 case (126). 
Multiplying both sides by and summing up over n from 3 to oo one gets 


■^T. 3 {s,t,z) = 2s‘^ [ dx I dyy{l-x) (^ 2 ( 5 , ^) + S 2 (f' 1 ^7 '^) \t'=tx-\-yu (134) 


0 JO 
00 


„1 00 

/ dx x^(i—x)^ — 

Jo 


1 dP 

z)\t>=-sx? (tsx(l - x))^’. 


pKp + 2)! dfP 


Using now that 

^ 3 ( 5 , t, z) = Si(s, t, zY - S' 2 (s, t)z‘^ + S'i(s, t)z, ^ 2 ( 5 , t, z) = Si(s, t, z) + Si(s, t)z, 

^S 3 (s,f,z) = ■^T,i{s,t,z) -282(3,t)z + Si{s,t), ‘282(3,t) = 25^^ (8i(s,t') + 8i(t',s)) 
one gets 




2 ;) =-^ + 2s^ I dx I dyy(l-x) (J:(s,t', z) + J:(t', s, z))\t'=tx+yu (135) 


„i 00 

/ dxx^(l —x)^y^ — 

Jo 


0 JO 
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p^pKp + 2)! dt'P 


dP 

( —(E(., t', z) + S(f', z))Y=.,xf (t3x(l - x)Y. 


And analogously for Zl=10 


»i 


^S(s,f, 2 ;) =-^ + I dx I dy y‘^(l-xf (J:(s,t\z)+ J:(t',s,z))\t'=tx+yu (136) 


dz 


„i 00 

/ dxx^(l —x)^y^ — 

Jo 


0 JO 
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p^pKp + 3)! dt'P 


dP 

( —(E(., t', z) + S(f', z))Y=-sxr (t3x(l - x))U 


The same equation with the replacement sot can be derived for S(t, s, z). 
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8 Conclusion 


Summarizing the obtained results one should admit that the calculation of the leading 
poles of the multiloop diagrams in nonrenormalizable theories is not a simple task con¬ 
trary to the renormalizable theories where they are given by the renormalization group 
equations. Nevertheless, we succeeded in writing down the recursive equations which al¬ 
low us to calculate the desired poles at any loop order. In a sense, these equations replace 
the RG ones for the nonrenormalizable case [22] • The difference is that in this case they 
are not algebraic but contain the integration over the Feynman parameters. 

Summation of perturbation series for the leading poles is also more complicated though 
the qualitative behaviour resembles that of renormalizable theories with the obvious re¬ 
placement —)■ g‘^s or g'^t. To get the sum of the inhnite series, we derived the differential 
equations and the task is reduced to the problem of hnding the hxed points and investigat¬ 
ing their stability properties. This is more explicit in D = 6 and looks more complicated 
in D = 8 ,10 due to the nonlinearity of the equations. 

In renormalizable theories the UV divergences either cancel in each order of PT like 
in the A/” = 4 D=A SYM theory, or are absorbed into the renormalization of the couplings 
and helds. In the nonrenormalizable case one may hope to get either the cancellation at 
each order (see the attempts in A/" = 8 SUGRA in 0=4) or the hnite limit as a result 
of summation of the inhnite PT series (like in D=4 QED where one has a zero charge 
behaviour). We have demonstrated that in higher dimensional maximally supersymmetric 
theories one may hope that the second opportunity is realized though the limit depends 
on the kinematics and we did not hnd the one where all the partial amplitudes are hnite. 
It might be that the D = 8 and H = 10 cases are better from this point of view. This is 
still to be found. 

The D = 6 SYM theory sometimes serves as a toy model for gravity since it has 
a dimensional coupling and the UV behavior similar to D = 4 gravity. The results of 
our analysis show that the hniteness of a theory which is our main goal is not reached 
so far. At the same time, it looks as if the loop by loop cancellation program does not 
work. It seems that the leading poles do not reveal any additional hidden symmetry, the 
dual conformal invariance being exploited already. Hopefully, the obtained new recursive 
relations and the differential equations for the inhnite sum of PT open the promising 
opportunity for the analysis of the UV divergences. 

9 Appendix A: The Master Integrals up to 5 loops 

All the diagrams up to 5-loops, which contain the leading pole, are presented in FigJT4| 
and FigjT^ The diagrams which can be obtained by the exchange sot are not shown. 
The rest of the diagrams can be found in [12]. The combinatorial factor is given on the 
left side of the diagram. 
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Figure 14: All the diagrams containing the leading pole from one to four loops 
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